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I. INTRODUCTION 

' 

' This work describes a set of lectures on Defect Structures in Field Theory. Although we focus mainly on issues 
, concerning scalar fields, we also deal a little with other topics. We work mainly in (1,1) dimensions, in flat space-time 
■ with metric such that diagiri^,^) = (1, —1), and we set the Planck constant h and the speed of light c to unit, and we 
use dimensionless fields and coordinates. 

The interested reader can find motivations and general investigations in the books [1-3], where the authors explore 
' several aspects of defect structures such as kinks, vortices and monopoles. There is another book, Ref. [4], in which 
. the reader may read important original articles in the subject until 1984. 
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where V{4i) is the potential, which specifies the model under consideration. 

We are working in (1,1) dimensions, and so the equation of motion which follows from the above model is 

If the field is constant we get dV/d(j) = 0. However, for static field configuration we have 



dx^ d(f> 

We are interested in solutions with finite energy. For the trivial constant solutions we have 
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Eo^ / dxVi^) (4) 



and for static solutions 
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In this last expression, the two terms identify the gradient and potential portions of the total energy, respectively. 

The trivial constant solution is represented by (j). It obeys 1^'(0) = and we should also impose that V{(f>) = 0, to 
make 4> to describe a null energy solution. On the other hand, static solutions must obey the boundary conditions 

lim ^^0 (6) 

x^-oc dx 
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and 

Urn (}>{x) 4> (7) 
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These conditions are necessary to make both the gradient and potential energies finite, independently. This behavior 
can be verified with the equation of motion: we rewrite Eq. (3) to get 

? = ±V2V (8) 
ax 

which also shows that the potential for a static solution cannot be negative. 

The behavior of the solution in the limit x oo may be equal to the former case, that is, ^(a; — » oo) (j); or it 
may be different, giving (^(a; — > oo) — > for (f being another trivial solution that obeys V'{(p) = 0, with V{(fi) = 0. 
These two possibilities can be described by the topological current 

Jt = le^'d.cl^ (9) 
which is conserved, that is, which obeys d/^jj, = 0. For static solutions the topological charge is given by 

Qt = ^(t>{x ^ oo) - ^(t){x -oo) (10) 

This charge sees the asymptotic bahavior of the static solutions, and it specifies two distinct cases: topological 
solutions, for which (^(x — > — oo) ^ (t>{x oo), and non-topological solutions, for which (j){x — oo) = <p{x oo). 
Thus, in general we can have topological or kink-like solutions, and non-topological or lump-like solutions. Wc can 
also have ribbons, which are kinks embedded in (2, 1) dimensions, and walls, which are kinks embedded in (3, 1) 
dimensions. 

We notice that there are several distinct ways to define the topological current, and this can be used to define the 
topological charge appropriately. 

In the case of a single real scalar field we can define the standard model, which engenders spontaneous symmetry 
breaking, with the potential 

VicP) = lil-ct>'f (11) 

This is the 0^ model, and it is plotted in Fig. 1. 
The equation of motion is 

and for static solutions we get 

= m^' - 1) (13) 

We are searching for extended but localized solutions, with finite energy. The equation above has two trivial solutions, 
given by <^± = ±1, which have zero energy and identify the classical vacua on the model. There are static solutions, 
given by 

(f)±{x) = ± ta,-nh{x) (14) 

They identify the kink (+) and anti-kink (— ), with center at the origin, a: = 0. In fact, the center of a static solution 
is arbitrary, because the model has translational invariance. 

The kink above has amplitude and width equal to unit. However, if we write the potential in the form 

V{^) = \xHa'-<pr (15) 



the kink would be (p{x) = atanh(Aaa;), and it would have amplitude a and width inversely proportional to Xa. 



3 



FIG. 1: Potential for the (A" model 



The energy for static solutions has the form 

E= I dxe(x) (16) 



where e{x) is the energy density, which can be written as 



We see that for the trivial solution <j), which satisfies V{(t)) = 0, the energy is zero. For the kink or anti-kink the energy 
is = 4/3, and the energy density is 

e(x) = sech^(a;) (18) 
We notice that the gradient and potential portions of energy are given by 

and 



Ep= dx V{(j)) (20) 



They vanish for the trivial solutions 0± = ±1. For the kink we get 

1 -4, 



ep(x) = eg{x) = -sech (x) (21) 



which shows that Eg ^ Ep. 

We can consider another model, described by the potential 

Vid^) = \cp'' ~ (22) 

which is plotted in Fig. 2. We see that this new potential is similar to the former one, but inverted; thus we name it 
the inverted t/)^ model. 

In this case, the equation of motion for static solutions is 

0.^-2*» (23, 
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FIG. 2: Potential for the inverted cj>'^ model. 

There is a single trivial solution, (f> — 0, for which ¥{(()) — 0. And there are static solutions, given by 

(/)±(x) = ±sech(a;) (24) 
They are lumps, and they have energy density in the form 

e{x) = sech^(x) tanh^ix) (25) 
In Fig. 3 we plot kink and lump, to compare their profile. 

<t>(x) ■ 




FIG. 3: Plots of the solutions for kinks and lumps. 

This last result shows that the energy density of the lump vanishes at the origin, the center of the lump, and is 
maximized at the points x = ±arctanh(\/2/2) « ±0.88. This suggests that lumps may have internal structure. In 
Fig. 4 we plot the energy density of both kink and lump, as given by Eq. (18) and (25). There we see that the energy 
density of the kink gets concentrated at the origin, but this is not the case for lumps. 

A. BPS States 

There is an alternative way to search for defect solutions. It is due to Bogomol'nyi [5] and Prasad and Sommerfield 
[6] . The point is that for non- negative potentials we can write 

VW = \wl (26) 




FIG. 4: Energy density for kinks and lumps, depicted as solid and dashed lines, respectively. 



where W = W{(j)) is a smooth function of the scalar field, and = dW/d(j). In this case the equation of motion for 
static solutions is 



It is possible to see that solutions of the first-order equations 



(27) 



(28) 



solve the second-order equation of motion. The solutions of the first-order equations are named BPS states, and it 
is possible to show that they are stable, minimum energy static solutions. In the present case, it is possible to show 
that the model only supports BPS states [7, 8]. 

To see how the BPS states appear, let us investigate the energy of the static solutions. We have 



2 



dx 



dx 



(29) 



We can write 



E — Ebps + 2 



(30) 



where 



Ebps - 1 1 W^[0(oo)] - W[c^{-^)] \ \ (31) 

is the bound in energy, which is attained for solutions that obey the first-order Eqs. (28). 

The BPS states present interesting properties. First, they are solutions of first-order equations, and their energies 
are easily obtained, since they only depend on the W{(j)), and on the asymptotic behavior of the field configuration. 
Second, the gradient and potential portions of their energies are equal, that is Eg = Ep = {1/2)Ebps- Another fact 
is that the BPS states appear very naturally in the supersymmetric environment, and there they partially preserve 
the supersymmetry. In the presence of supersymmetry, W is named superpotential; see Refs. [9-14] for more details 
on the subject. 

There are different motivations to investigate models which support BPS states. In particular, in Ref. [15] one 
investigates elastic deformations in liquid crystals, and there the decomposition in terms of BPS states has led to 
wide class of solutions. Another interesting work was recently done in supersymmetric quantum mechanics [16]. In 
this case, it was shown how to understand shape invariance in terms of the BPS structure of the system. 
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B. Linear Stability 

Linear or classical stability of static solutions are investigated as follows: for ^ = (p{x) being static solution, we can 
write 

(l){x,t) = (l){x) + r]{x,t) (32) 

where r]{x, t) describes fluctuations around the solution (l){x). We suppose that the fluctuations are small, and we only 
consider linear contributions in r/. We then use the above conflguration and the equation of motion (2) to get, up to 
first-order in rj, 



d^T] d'^r] (PV 



r? = (33) 



Since the static solution only depend on x, we can consider 

r]{x, t) = ^ r]n{x) cos{wnt) (34) 



to obtain the Schrodinger-like equation 



where 



H rjn {x) = wl r]n {x) (35) 



d2 d^V 

H = — — . + 



(36) 



In the Schrodinger-like equation, we notice that the presence of negative eigenvalue transforms the cos factor in 
(34) into a cosh, and this would violate the supposition that we are dealing with small fluctuations. Thus, we need 
that the associated quantum-mechanical problem have no negative eigenvalue, to make the solution stable. 

An interesting fact is that the associated quantum-mechanical problem always contains a zero-mode, that is, an 
eigenfunction with zero eigenvalue. We prove this assertion by noting that if (f){x) is a static solution, we use the 
equation of motion and (j)' = d(j)/dx to obtain 

which is a Schrodinger-like equation with zero eigenvalue. Thus, we identify the zero-mode with the derivative of the 
classical solution itself. The presence of the zero mode is a general fact, which appear in decorrence of the translational 
invariance of the model under investigation [17] . For instance, for the (f)^ model the kinks are (j)± [x) = ± tanh(x) . They 
are stable, and the zero-modes arc </''-j-(a;) = ±scch^(a;), as wc are going to show below. For the inverted c/)'^ model, the 
lumps are 0±(a;) = ±sech(x). They are unstable, with zero- modes as 4'±{x) = =F tanh(a;) sech(x). Their instabilities 
appear from the fact that such zero-modes cross zero at the origin, so they have nodes and cannot represent lowest 
bound states [18]. 

In the case of the non- negative potentials described in (26), there are BPS states, which obey first-order equations. 
They are all stable, and this can be shown as follows: the Hamiltonian for BPS states are given by 

H = -^+W}^ + W^ (38) 



and now we can introduce the first-order operators 



and they allow introducing the Hamiltonians 



-l±W,, (39) 



H± = SiS± 



(40) 
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which can be written in the general form 

H± = -^+ Wl^, ± W^W^^^ (41) 

These Hamiltonians are supersymmetric partners [19, 20] in quantum mechanics. Now, since Hj^ is the Hamiltonian 
of equation (38), we see that it is non- negative, thus it has no negative eigenvalue, and this ensure hnear stabihty for 
the BPS states. We can also see that the zero- mode r}o{x) obeys 

^ = ±^..%(x) (42) 

and it can be written as 

T1o{x)=Ae^S'^'=^i'f (43) 

where (j) = (t^{x) stands for the static solution. The constant A is used to normalize the wave function, but this only 
occurs for one of the two signs of the exponential. There is another way to write the zero-mode; it is 

r?o = (44) 

where <j) = (j){x) is the static solution, the BPS state. The zero-mode ?7o(a;) is also known as the bosonic zero-mode, 
since the inclusion of fermions in the model will also lead to fermionic zero- modes [21]. 

We consider the cf)^ model. The topological solutions are (f)± (x) = ± tanh(x) and the potential for the corresponding 
quantum-mechanical problem has the form 



U{x) = 



cPV 



= 4-6sech^(a;) (45) 



This is the modified Poschl- Teller potential, which was investigated long ago [22, 23]. We consider the potential 

U{x) = a + b tanh(a;) — c sech^(x) (46) 

with a, b, c real and positive parameters. For a > b and a < c, the continuum spectrum presents reflecting states, for 
energies in between a — b e a + b, and other states for energies higher than a + b. If 6 is zero, the continuum does not 
contain reflecting states. For a > b and a < c, the discrete spectrum may contain the zero-mode and other states, 
with positive and negative eigenvalues. In fact, the bound states present eigenvalues given by 



6^ 

4^2 



(47) 



where 



and n = 0, 1, is such that 



^n = \/c+i-n-i (48) 



n<Vc+i-i-y- (49) 

In Ref. [22] there is a misprint: the last term of the result for the bound states there obtained has to be multiplied 

by the factor 1/4. 

The above results can be used to show that for the c/)'^ model, the potential (45) leads to two bound states, one being 
the zero mode, and the other having eigenvalue 3. For the inverted c/)^ model, the potential for the quantum-mechanical 
problem is 

?7(x) 1 - 6scch2(x) (50) 

It has the very same form of the potential for the (j)'^ model, but it plots differently, in a way such that its zero mode 
corresponds to the higher bound state, thus introducing a negative eigenvalue, —3, which makes the lump unstable. 
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C. Examples 



We now consider some examples, to illustrate the above results. The first example is the sine-Gordon model, which 
is described by the potential 



Here we can write W{(j)) = — cos(^). The first-order equations are 



dx 



= ± sm{(j)) 



(51) 



(52) 



The minima of the potential are given by jio = 0,_^±i = ±7r, etc. There is an infinity of topological sectors. They are 
all equivalent, and in the sectors 4>-i ^ 0o and 4>o ^ the solutions can be written as 



01 (x) = ±2arctan(e±^) 

The energy is E — Ebps — 2. 

Stability investigations lead us to the quantum mechanical potential 

Z7(a;) = l-2sech^(a;) 

which has just one bound state, the zero mode, with normalized eigenfunction given by 



mix) 



- sech(a;) 



(53) 



(54) 



(55) 



There are other sectors, and there are other solutions; see, e.g., Ref. [1, 4]. 

Another example is the double sine-Gordon model. It may be described by the potential 



v{4>) 



-(r + cos((/)))2 



(56) 



where r is a real parameter. This model was recently investigated in [24]. It can be described by the superpotential 



W(<A) = 



--{rcj) + siTL{(j))) 



(57) 



We consider r e (0, 1) and the singular points, dW/d(f> = 0, identify the minima of the potential. They are periodic 
and, in the interval — 27r < (j) < 27r, we can identify cj) = ±7r ± arccos(r). We see that the maxima in ±7r and the 
minima in ±7r ± arccos(r) collapse to two minima in ±7r in the limit r ^ 1. This fact is illustrated in Fig. 5, where 
we plot the double sine-Gordon potential for the values r = 1/3, 2/3 and r = 1. 

This model contains two topological sectors, one large and the other small, with energies 



E' 



BPS 



4\/l - r + 4r 



and 



E' 



BPS 



4vr 



TT — arccos(r) 
VT+r 

arccos(r) 



4r 



(58) 



(59) 



We notice that E'-^pg 



Egpg + Anr/y^l + r, and in the limit r ^ 1 we get Ep^g 
Stable solutions which solve the first-order equations are given by 



2\/27r e Ei 



BPS 



0. 



and 



(x) = ±2 arctan 



. {x) = ±7r — 2 arctan 



1 + r 



tanh(-\/l — rx) 



1 + r 
1-r 



tanh(\/l — rx) 



(60) 



(61) 
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FIG. 5: Potential for the double sine-Gordon model, for r — 1/3, r = 2/3, and r — 1, depicted with dot-dashed, dashed and 
solid lines, respectively. 



It is interesting to see that the double sine-Gordon model can be used to map specific magnetic materials, like the one 
studied in Ref. [25]; there one identifies three distinct phases: anti- ferromagnetic, canted, and weakly ferromagnetic. 
These phases can be described by the above model, in the cases r = 0, 0<r<l and r = 1, respectively. Moreover, 
in the canted phase, there also appear two distinct defect structures, the large and small solutions. 

Another example is described by the (p^ model. Here the potential has the form 

We have W{<j)) — (1/2)0^ — (l/4)(/)^, and the minima are 4>o — and 4>± — ±1. There are two topological sectors, 
equivalent. The solutions are 



(/>|(a;) =±^i[l±tanh(a:)] (63) 

and they identify the kinks and anti-kinks in the two sectors. Their energies are given hy E = Ebps = 1/4- 
Stability leads us to the quantum-mechanical potentials 

5 3 15 

U±{x) ^ -T ^ tanh(x) - — sech^(a;) (64) 

The terms which contain the hyperbolic tangent make the potentials asymmetric. This introduces reflection states 
into the continuum, although there is only one bound state, the zero-mode; see Ref. [26] for more details. 
Another model is given by [27] 

V{4>) = \{l-m' (65) 

It is similar to the 4)^ model. It has minima in (}>± = ±1, and we can write W{(l)) — (f> — (1/2)|(/)|(/). The first-order 
equations are 

^ = ±(1-101) (66) 



The solutions are 



, , tanh(a;/2) , , 

' 1 + tanh .x /2) ^ ^ 



Their energies are given hy E = Egpg = 1. In Figs. 6 and 7 we plot potentials and kinks, for the above model, and 
for the 0* model, for comparison. 
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FIG. 6: The potential of Eq. (65), depicted with dashed line, and the potential of the model. 




FIG. 7: The kink of Eq. (67), depicted with dashed line, and the kink of the 0* model. 
Stability investigations lead us to the potential 

Uix) = 1 - 25{x) 



(68) 



which is attractive. There is just one bound state, the zero mode. To see this explicitly, note that the Hamiltonian 
can be written as 



Thus, we see that the zero-mode obeys 



+ 1 -2 Six) 



d \x\ 



d \x\\ f d \x\ 
dx X I \ dx X 



rioix) = 



This equation can be solved easily, given the normalized zero-mode 

r]a{x) = exp(-|x|) 



(69) 



(70) 



(71) 



The above model was investigated in [27, 28], leading to exact solutions for a model described by coupled linear 
chains. More recently, it was studied in [29, 30] with other motivations. In particular, in [30] we have introduced 
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FIG. 8: The potential of the Eq. (72), depicted with dashed line, and the potential of the (f)^ model. 




FIG. 9: The kink of the Eq. (73), depicted with dashed line, and the kink of the (j>^ model. 

another model, inspired in the above model and in the (j)^ model. It is described by 

T/(^) = i<^2(l- 101)2 



(72) 



There are minima in cf) — and (j)± = ±1, and we have W{(l)) = (l/2)(?!)2 - {l/3)\(j)\(j)'^. There are two topological 
sectors, equivalent, with energies given hy E — Ebps = 1/6. The solutions are 



b^ix) = ±-[1 + tanh(a;/2)] 
b^ix) =±i[l-tanh(a;/2)] 



(73) 
(74) 



They are similar to the solutions of the (p^ model; see Figs. 8 and 9 for a comparison. 

An interesting feature of this new model is that the topological sectors are symmetric, related to the local maxima 
they have. This makes the topological solutions uniform, leading to a symmetric quantum-mechanical potential which 
contains no reflection states, as it appeared in the 0® model. This model is of direct interest to study first-order phase 
transition, and we believe that it may be of some use in applications to Condensed Matter. 

We investigate stability of the solutions of the above model. We get to the potential 



Uix) 



1- -sech2(a;/2) 



(75) 



12 



The Hamiltonian can be factorized as 

f- 

dx J \dx 



+tanh(a;/2)) f + tanh(x/2) ) (76) 



Thus, the zero-mode should obey 



^ + tanh(a;/2)^ 770(0;) = (77) 
We integrate this equation to obtain the normahzed zero-mode 

= \^sec\i\x/2) (78) 

There is another bound state, with eigenvalue 3/4, and no reflection state appears in the continuum. 
Another example is given by the potential [31] 

y(<^) = \^cc\,\<i>) (79) 

The novelty here is that the potential has no minimum, but it supports topological defects, as shown in Refs. [31, 32]. 
Our participation in this problem started with [33], where we have shown that this model engenders first-order 
equations given in terms of the superpotential 

W{(j)) = arctan[sinh(0)] (80) 

The first-order equations are 

^=±sech(<^) (81) 
dx 

They are solved by 

(p±{x) = ±arcsinh(a;) (82) 
These defects are different: they are very diffuse, as we see from the energy density 

ei^) = ^ (83) 

However, the energy is finite, given by = Ebps = 

It is interesting to see that if we write the topological current in the usual form j^^ = {\/2)e^^d^(l>, we would have 
problem to calculate the topological charge. However, since the topological current is not unique, we can write 

ji" = e'^'d^Wicp) (84) 

In this case, we are relating the topological charge with the energy of the solution, which seems to be appropriate, 
since the energy should be well-defined. 

We investigate stability of the above solutions to get to 

This is a volcano-like potential, as we show in Fig. 10. It has a single bound state, the zero-mode, with normalized 
eigenfunction given by 



In this case, there is no gap between the zero-mode and the continuum. 
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U(x) 




FIG. 10: The volcano-like potential. 



Another model is described by the potential [34] 



(87) 



The novelty here is that this potential is unbounded from below, and it has two inflection points, at the values (j) = ±1, 
such that V{<j)) = 0. However, there are topological solutions connecting these inflection points. They are given by 



(88) 



The appearance of topological solutions connecting inflection points was first shown in [34]. 
We investigate stability to get to the potential 



U{x) = 12- 



-l/4 + a;2 
(l + a;2)2 



(89) 



which is also volcano-like. It has only one bound state, the zero-mode, and there is no gap between the zero-mode 
and the continuum. The zero-mode is given by 



1 



37rY (l + x2)3 



(90) 



In the last two examples, there is no gap between the zero-mode and the continuum, and so we should investigate 
how the quantum corrections should appear, to see if there are subtleties, as it occurs in models where a gap is present; 
see, e.g., Refs [12-14] for more information on this issue. 

Another model is described by the potential [35] 



V{4>) = i(0(p-i)/p-0(p+i)/p)2 



(91) 



where p is a real parameter. For simplicity, here we consider p odd integer, p = 1, 3, 5, ... This model can be seen as a 
generalization of the model, which is obtained for p = 1 . 
In this case we can write 



WU) = ^^0(2f-l)/P ^^(2P+1)/P 

2p — 1 2p + 1 



and the first-order equations are solved by 



4>±{x) = ± tanh''(a;/p) 



(92) 



(93) 
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FIG. 11: The 2-kink solution for p = 3 and 5, depicted with sohd and dashed hnes, respectively. 

The energies are given by E^pg — 4p/{Ap'^ ~ !)• In Fig- H we plot some solutions, and there we see that they connect 
the minima ±1, passing through the minimum at zero with vanishing derivative. 

As far as we can see, this is the first time that topological solutions connect non-adjacent minima, passing through 
an intermediate minimum. This fact is only possible because the minimum at zero is not a true minimum, that is, 
it is not good for quantization because it cannot be used as a perturbative ground state, since the mass squared at 

= is divergent. For this reason, the above solutions are named 2-kink solutions. 

The above solutions present normalized zero-mode, described by 



= J^^^ tanh(P-i) (^^ j sech^ (^^ ] (94) 



They have maxima at the points 



/n - 1 

x± = zbparctanhy — j-y (95) 
which can be used to identify how the two kinks localize inside the 2-kink solution. 

III. TWO-FIELD MODELS 

Let us now consider models described by two real scalar fields, given by 

£ = ia^09^0 + ia^x9^x-mx) (96) 

where V{(l),x) is the potential, which specifies the model to be studied. In (1, 1) dimensions, the equations of motion 
are 

fi2j, pi^r 

(97) 
(98) 

(99) 
(100) 



and for static solutions we have 



5^0 




dV 


dp 










dV 






dx 






dV 






90 




d\ 


dV 




dx^ 


dx 
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Non-linearity makes it hard to deal with these equations on general grounds. Some investigations were done in 
Refs. [36-38]. In particular, in [36] one illustrates the difhculties to solve Eqs. (99) and (100). This problem was 
further considered in [39] and there it was proposed a method, named trial orbit method, which could perhaps be of 
some use to solve the above equations. In 1995 we started our investigations on models described by two real scalar 
fields. We started with the potential [40] 



In this case, the equations of motion for static fields can be written in the form 

dx 



y{M = 7^'^l^7^K (101) 



^-W^W^^^W^W^^ (102) 



^^^-W^W^^ + W^W^^ (103) 



and the energy associated with these configurations could be written as 

2 / , X 2 



This expression can be rewritten in the form 

2 



(104) 



(105) 



where we have set 

Eb = Ebps = \AW\ (106) 

with 

AW = iy((/)(oo), x(oo)) - W{(j){-oo), xi-oo)) (107) 

This procedure shows that the energy is minimized to i? = Ebps, for field configurations which obey the first-order 
equations 

g=±H^, (108) 



dx 
dx 



On the other hand, if we differentiate these equations with respect to x, we can see that their solutions also solve the 
second-order equations of motion for W^-^ = W-^^, which is the condition for W to present second derivative, as the 
potential (101) asks implicitly. 

It is interesting to see that solutions of the first-order equations have the same portion of gradient and potential 
energies. The first-order equations can be seen as a dynamical system, and so wc can use all the mathematical tools 
available for dynamical systems to deal with these equations. In particular, the algebraic equations Wrf, = and 
= can be solved to give the set of singular points of the respective dynamical system. In the field theoretical 
model, this set of singular points constitute the absolute minima of the model. We can characterize the absolute 
minima with Vi = {(pi, Xi), ^ = 1)2, ... Thus, each pair of minima Vi e Vj may represent a topological sector, which may 
be BPS or non-BPS. 

A given sector ij is BPS if the energy is not zero, that is, if Wi = W{(f)i,Xi) is distinct from Wj = W{(f)j, Xj), 

to make the energy -E^ps = l^^ijl — ~ non-vanishing. On the other hand, if Wj = W(0i,Xi) is equal 
to Wj = W{(l)j,Xj), there is no BPS energy available, and this shows that there is no BPS solution connecting the 

minima that identify the sector. In this case, however, the sector may have configurations which solve the equations 
of motion. These sectors are usually named non-BPS, because the field configurations do not solve the first-order 
equations. 

In the BPS limit, models described by two real scalar fields are simpler, because we have to deal with pairs of 
first-order equations. We then focus on the equations (108) and (109). We can write 

Wxd(l)-W^dx = (110) 
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This equation is easily solved when W{4>, x) is harmonic. In this case, the solutions describe an orbit in the (0, x) 
plane, given by F{(f),x) = 0, which is obtained with 

OF dF , , 

- = W,, - = -W, (111) 

In Refs. [7, 8] we have shown that when W is harmonic, all the solutions are BPS solutions. 

When W{(f),x) is not harmonic, we have to search for the integrating factor 7(0,%). This factor is needed to make 
another function, F{(f),x), harmonic, such that 

^=I{cj,,x)W^, ^ = -I{^,x)W^ (112) 

The main problem with the integrating factor is that we do not know how to get it in general, and so the trial orbit 
method is yet of good use to find solutions in specific models. 

The trial orbit method is very effective when it is applied to first-order equations. This was shown in Ref. [41], and 
the main steps can be summarized as follows: 

First step. Choose a BPS sector. This is done with two minima, Vi = {4>i,Xi) a^^d = {'t'j 1X3)1 such that 

Wi = W{<j>uXi)^W^ = W{<l>j,Xj)- 

Second step. Choose an orbit. This is done with a given function Fij {(f>, x) = 0, which is compatible with the 
minima that specify the sector, which satisfies Fij{(pi, Xi) = and Fij{(j)j, Xj) = 0. 

Third step. Try the chosen orbit. This is done taking the derivative of the chosen orbit to get 

^H', + ^Vr,=0 (113) 

If this new expression is compatible with the chosen orbit, it is a good orbit and can be used to decouple the first-order 
equations. 

These three steps are very efficient to solve the first-order equations, to find BPS solutions, as we illustrate below. 



A. Linear Stability 

We can investigate linear stability of the BPS solutions in a general way. We write </> = 0(a;, t) and x = xi^i '^s 

<l){x,t) = (f){x) + ij{x,t) (114) 
X{x,t) = x{x)+^{x,t) (115) 

where rj{x,t) and £,{x,t) are fluctuations around the static configurations (j)(x) and x{x)i which solve the first-order 
equations. Since the classical solutions are static, we can write 

?g::0=i;(ii1)™("'-" 

We use these expressions into the equations of motion (97) and (98) to get to the Schrodinger-like equation H^n{x) = 
wl^^nix), where 



and 

We consider the potential (101) to get 



*n(a;)=fHS ) (118) 



= W^^ + W^W^^^ + W^^ + W^W^^^ (119) 
Vxx = + W^W^xx + Kx + WxWxxx (120) 
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and 

V'^x = yx4> = W^^W^^ + W^W^^^ + Wl^W^^ + W^W^^^ (121) 
In Ref. [42] we have shown how to introduce first-order operators such that 

which lead us with the Hamiltonians 

H± = si S± (123) 

Now, in (117) we see that H is just the partner of The Hamiltonians e H- can be seen as supersymmetric 
partners in quantum mechanics. The factorization of H impUes that it is non-negative, and this proves stabiUty of 
the BPS solutions in general. Moreover, the zero-mode solves the equation 



dx 1^x0 W^x 



^0 (124) 



and so it has the form 



where A is the normalization constant. This result can be seen as a natural extension of the former result, for a single 
field, given by Eq. (44). 

B. Example 

We illustrate the above investigations with an example. We consider the function 

m</',x) = 0-^<^'-r-</'x' (126) 
where r is a real parameter. The potential has the form 

VicP, x) = l{l-cl>^?-rx^+r{l + 2r)<p\^ + \r\^ (127) 

This model was first studied in Ref. [40]. See also Refs. [11, 43-52] for related investigations. 

The potential (127) is plotted in Fig. 12. For static solutions the equations of motion are given by 



The first-order equations are much simpler 



-2(/) + 2r(l + 2r)(l)x^ + 2<p^ (128) 
-2rx + 2r(l + 2r)(/>^x + 2r^x^ (129) 



^ = ±1^^2^rx^ (130) 

^ = T2r^X (131) 

The minima of the potential are (±1,0) for r < 0, or (±1.0) and (0, ±l/^/r) for r > 0. We consider the richer case, 
for r > 0. Hero there are four minima, given by vi = (1,0), V2 = (—1,0), V3 ~ (0, l/-\/r) and 114 = {0, —1/^/r). There 
are six topological sectors: sector 1, which connects the minima vi and V2; sector 2, which connect the minima W3 
and V4; sector 3, which connect the minima Vi and V3; sector 4, which connect the minima Vi and V4; sector 5, which 
connect the minima V2 and V3; and, finally, sector 6, which connect the minima V2 and V4. 
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The value of W at the minima are given by: Wi = W{1, 0) = 2/3, W2 = W(-l, 0) = -2/3, W3 = W(0, 1/V^) = 0, 
and W3 = W(0, — l/Vr) = 0. Thus, there are five sectors which are BPS, one with energy E^pg = 4/3 and the other 
four with the same energy E^pg = Ej^pg = E^pg = E^pg = 2/3. There is one sector which is non-BPS; it has 
energy E^^pg, which can only be obtained when one shows the corresponding topological solution explicitly. 

To find explicit solutions, let us first consider solutions which depend only on a single field. We first set x = in 
the equations of motion. We get 

0=-2</,(l-</,2) (132) 

which is the equation for the 4>* model, which has the solutions 4'±{x) = ±tanh(a;), with energy 4/3. Thus, the pair 
of solutions 

(?!)±(x) = ±tanh(a;), x = (133) 

are solutions for the above two-field model. In field space, these solutions describe a straight line connecting the 
minima vi and V2, which identify the topological sector 1. We also see that the imposition x = reduces the 
first-order equations to the equation 

f=±(I-^^) (134) 

which have the same solutions, showing that they are BPS solutions. 
On the other hand, if we set = in the equations of motion, we get to 



which is solved by 



"''^ =-2r2x(l/r-x') (135) 



X±(a;) = i-^ tanh(v^a;) (136) 



with energy E^^pg = 4/3\/r, which depend on r. Thus, the pair of solutions 



' = 0, x±{x) = tanh(^a;) (137) 
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represent straight line orbit which connects the minima and v^. In this case, however, the solutions are non-BPS, 
because they do not satisfy the first-order equations. We notice that the energy of the non-BPS solutions depend on 
r; for r G (0, 1) it is higher than the energy of the BPS states in sector 1. 

The above model has an interesting feature. For r = 1, we can rotate the {4>,x) plane to 0i = (j^ + x and (j)2 = (f> — Xj 
to see that (f>i and 02 decouple, describing two single field models, instead of a single model of two fields. See Ref. [52] 
for other details on the issue. 

In order to illustrate the trial orbit method, let us investigate solutions in the several BPS sectors of the model. 
First, we consider sector 1, identified by the minima vi = (1,0) and V2 = (—1,0). We consider an elliptic orbit, 
described by a<^^ -I- bx^ = c, with a, 6, c being real parameters which identify the orbit. This orbit must contain the 
two minima which identify the sector 1; these conditions eliminates two parameters, and we get to the more specific 
orbit 4>^ -\- bx^ = 1. We now derivate this orbit and use the first-order equations to get to t^^ + r(l + 2b)x^ = 1, which 
is compatible with our orbit for b = r/{l — 2r). Thus, the orbit 



1 



■X 



1 



(138) 



is a good orbit for r G (0, 1/2), and can be used to solve the problem. We use this orbit to write 



We use this into the first-order equation for to obtain 



which are solved by 

Thus, for the other field, Xj we get 



-{x) — ± tanh(2ra;) 



X±(x) = ±\/ 2 sech(2r2;) 



(139) 

(140) 
(141) 

(142) 



There are four solutions for the pair (0, x)? and in Fig. 13 we plot a pair of solutions for r — 1/6. It is interesting 
to notice that the above solutions can be used to describe solitons in ferroelectric materials, as we have pointed in 
Ref. [53]. These solutions are BPS states. Their energies have the value 4/3, and they change to the solutions that 
describe the straight line orbit in the limit r 1/2, as it should be expected. 




FIG. 13: A pair of solutions for the two-field model. The solid and dashed lines describe the <j> and x fields, respectively. 

There are other orbits, in the other sectors, and they can be obtained as we did above. Orbits in the other sectors 
were first obtained in Ref. [46]. Also, an interesting investigation done in Ref. [47] obtained the integrating factor for 
the pair of first-order equations, solving completely the problem of finding BPS solutions for the model. 
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We can study stability of the solutions explicitly. We illustrate this point investigating the simplest case, concerning 
stability of the solutions for the straight line orbit in sector 1. In this case, the fluctuations in the fields (j) and x obey 
the equations Hr]n{x) — 77,1 (a;) and HS^n{x) = w'^^n{x), respectively. The Hamiltonians are given by 

H = -— +4-6sech2(x) (143) 

and 

^2 

H = — — +4r2-2r(l + 2r)sech^(a;) (144) 

The first problem is the modified Poschl- Teller problem. It supports two bound states, the zero-mode Wo = and 
wi = -y/S- The second problem is also of the modified Poschl- Teller type, but now the number of bound states depends 
on r, the parameter that controls the form of the potential. We have 

Wn = (4r — n)n (145) 

where n = 0, 1, 2, < 2r. Thus, if r e (0, 1/2], only the zero-mode appears. For r E (1/2, 1], there are two bound 
states, the zero- mode and another one, with eigenvalue wi — Ar — 1. We illustrate this in Fig. 14, in which we plot 
the potential of Eq. (144) for the values r = 1/2, r = 1 e r = 3/2. We notice that for x = 0, the investigations for 
solutions reduce to that of the (f)^ model; however, stability is very different, and this will certainly show up at the 
quantum level. 
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FIG. 14: Potentials for r = 1/2, 1, e 3/2, depicted with solid, dashed and dot-dashed lines, respectively. 



C. Junction of defects 

In two-field models, in general the potential depends on the two fields, and so the minima are in the {(l),x) plane. 
Thus, we can find situations where the minima constitute a non-collinear set of points, which is not possible for 
models described by a single field. For this reason, two-field models may give rise to junctions of defects, which are 
regions where two or more defects intersect, and this lead us to more general situations that describe junctions of 
defects, and the tiling of the plane with a network of defects. Some investigations on this subject can be found in 
Refs. [50, 51, 54-58] and also in the more recent Refs. [7, 8]. 

An interesting way to describe junctions of defects is with a complex field. We consider models described by the 
Lagrange density 

C^d^^d^^-VM) (146) 

with potential 



vM) = \w'm'' 



(147) 
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given in terms of an holomorphic function W = W{ip), with W'{ip) = dW/d<fi. For this model, the equation of motion 
for static solutions has the form 



(148) 



The minima of the potential are given by W'{ip) = 0. We suppose that there is a finite number of minima, given by 
Vk, k = 1,2, N, where N is an integer. We are interested in solutions which obey the boundary conditions 



ip{x -OO) Vi, 

where Vi is one in the set of N minima. 

We can show that solutions of the first-order equation 

dx 



dip 



-oo) 



(149) 



(150) 



also obey the equation of motion, for ^ being real and constant. 

The above model is special, since the equation of motion is equivalent to the set of first-order equations parametrized 
by the phase ^. To show this we use the first-order equation to write 



dip 
dx 



W'{ip) 



(151) 



This shows that solutions of the first-order equations solve the equation of motion. On the other hand, we can 
introduce the ratio 



R{ip) 



1 



dip 
W'{ip) d^ 



We use the equation of motion to get 

dR 
dx 



1 



W'{ip) 
1 



\W'{^)\' 



dip 



dx 



W"{^) 



(152) 



(153) 



We can also write 



and 



dp 



dip 



dx 



dx 



d_ 
dx 



dip 



dx 



W'{ip)W"{^)^ + W'{ip)\W"{^) 



dx 



dip 
dx 



The last two equalities lead us to 



(154) 



(155) 



S{ip) = \W'{^)\ 



dip 



dx 



(156) 



which does not depend on x, that is, which obeys dS{p)/dx = p solves the equation of motion. However, in the 
limit X — *■ —00 the field should become ip Vi, and dip/dx — + and W'{ip) 0. Thus, S{ip) is shown to vanish, and 
this implies that R{ip) is constant. Besides, since S{ip) vanishes, we can obtain = |iy((^)|. Thus, we can write 

|i?((^)[ = 1. and so R(p) = e'^, with ^ being real and constant. The last result is nothing more than an alternative 
way to write the first-order equation. This is general, and it shows that in this model all the static solutions are BPS 
solutions. A direct consequence of this result is that for models described by two real scalar fields, if = W{(t), x) is 
harmonic, then the model cannot support non-topological solutions. 
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We can consider the model [7, 8] 

= (157) 
where n = 0, 1, 2, and N = 2,3, ... are integers which characterize the models. The potential has the form 

y = ^(w)"(i-<^'^)(i-<p'^) (158) 

The models engender discrete symmetry, Z^, for every n, and they can be used to generate defect junctions, of 
intrinsic interest, and for tiling the plane with regular polygons. 
The case n = is special. It has 

WNiv) = ^-Y^v'^'' (159) 

which were already studied in [59], with other motivations. 

The general model has minima at vj^ = cxp(?;^^), for = 2Tr{k/N), k = 1,2, N, independently of n, and also in 
f = for n ^ 0. There are N minima in the unit circle, and the origin is also a minimum for n ^ 0, with multiplicity 
n. For the topological sectors, for n ^ there appear N radial sectors. The energies of the radial sectors are given by 

N 

^ (l + n){l+n + N) 
for n = 1, 2, ... and = 2, 3, ... The energies of the other sectors are 

E^^^ = 2E^^o\Mi^+n){kn/N)]\ (161) 

for n = 0, 1, ... and = 2,3, ... 

We notice that if the pair n and N obeys n = N —1, the energies E"^ ^ vanish for every A; = 1, 2, A^. In this case 

the model can only support radial sectors, with energies E^~^ = 1/2N. 

As one knows, to tile the plane with regular polygons [60, 61], we only have three possibilities: using equilateral 
triangles, squares, or regular hexagons. Now, if we ask for the most efficient way of tiling the plane, we see that 
it is given by the hexagonal tiling. Efficiency here implies the greater relation area/perimeter for the three regular 
polygonal tilings. Since the use of hexagons requires the presence of triple junctions, we then focus attention on the 
above model, for the case n = and A^ = 3. This case can also be of interest to nanotubes and fuUerenes, as explored 
in [62, 63]; see also Refs. [64-67]. 

We use n = and A'' = 3 to get 

W{^)=V>-\^'' (162) 

The three minima are Vi = 1, and 112,3 = 1/2 ± i-\/3/2. The energy of the defects are given by E^ = 3v^/4, and this 
should be the energy density of the junction in the regions far away from the junction. However, if we neglect the energy 
density related to the junction itself, we can write the total energy of the triple junction as E = 3LE^ = 9\/3L/4, 
approximately, for L being the size of each junction leg. 



IV. DEFORMED DEFECTS 



In the case of models described by a single field, recently we have been able [34, 68] to introduce a procedure, which 
we call deformation procedure, where we start with a given model, that support defect structure's, and we deform it, 
to write a new model, which also support defects that can be written in terms of the defects of the given, starting 
model. To make this idea practical, let us start with a model defined by the potential V{(f>), which support static 
solutions, given by 0(x). To deform this model we get a function / = f{(l)), which is well-behaved, invertible, which 
we name deformation function, and we write another potential, as follows 



(163) 
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We have two distinct models, given by 

£ = ^d^<f>d^^cl>-V{ct>) (164) 
C = ld^^d^^-V{v) (165) 

However, if V and V are related by Eq. (163), and if (p±{x) are defects of the first model, then the new model has 
defect solutions given by 

<fi±{x) = f-\<f>±{x)) (166) 

written in terms of the inverse of the deformation function. 

To prove the above result, let us consider the equation of motion for static solutions in the two models 

We require that the energies are finite, and that the solutions obey boundary conditions as before, in Sec. II. Thus, 
we can write 

In the first case, the energy density has the form 

2 / , , \ 2 



In the other case, we get 



1 / d(f\ ^ ~^ f dip 

2 [d 

We can use the last expression to relate both densities 

(d//#) 



2 



On the other hand, we can use (p = f{f) to obtain 



d'^cj) _ d / df dip 
dx^ dx \ dip dx 



dip dx"^ dip'^ \dx J 
Thus, we use d^ip/dx^ and dip/dx from the equation of motion to get 



dj d^ip d^f fdipY 

■±77^ + :7^ ±) (172) 



dx'^ dp dp dpP' 
However, since V{ipi) = V {(p) / {df / dip)"^ , we can write 



d'^cj) df dV ^-d'^f 
Y .1 ~ 2V (173) 



dV _ 1 rf^f 1 dV 

dip ~ ^"^^ {df/dp^f dip^ ^ (df/dip) d(P ^ ' 
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and so we get 



dV 



(175) 



which is the equation of motion for the first model. This resuh shows that sohitions of the second model are 
related to solutions of the first model. To complete the proof, we need to show the contrary, but this is similar and 
straightforward. 

It is interesting to realize that we can deform the deformed model once again, etc, and this procedure lead us with 
an infinity of models, as we illustrate in the diagram below. 




We notice that the deformation procedure works for models that support topological, kink-like, or non-topological, 
lump-like, solutions. However, for models which support BPS states, the deformation procedure is simpler, since now 
we can investigate the first-order equations. We have 



dx 



dW_ 
d(t) 



and 



The last equations can be written as 



d(p 
dx 



± 



1 dW 
df/dip df 



dl^^dW 
dx df 



(176) 



(177) 



(178) 



and now they are like the equations for the first model, with the change (j) f. Thus, their solutions obey f{f±) = 4>±, 
and so we get ip±{x) = f^^[4'±{x)\, in accordance with Eq. (166). 

We illustrate the procedure with an example. We take the standard 0* model, and the function /((^) = sinh(0). In 
this case the potential 



2\2 



changes to 



V{<P) = -{l-<i>') 



1 



V{(j)) = ^sech\(l))[l - sinh2(^)]2 
It is not hard to verify explicitly that the static solutions of this new model are given by 

4>±{x) = ±arcsinh[tanh(a;)] 
Besides, the energy density of these solutions are 

e{x) = cosh^[tanh(a;)]sech*(a;) 



(179) 



(180) 



(181) 



(182) 



where s{x) = sech^(a;) is the energy density for the 0* model. In Figs. 15 and 16 we plot defects and energy densities 
in both models, for comparison. 

Another example can be given; we consider the (j)'^ model, but now we use the function 



i/p 



(183) 
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FIG. 15: Defect solutions for the 0* model, and for the deformed model, depicted with solid and dashed lines, respectively. 




FIG. 16: Energy densities for the model, and for the deformed model, depicted with solid and dashed lines, respectively. 

where p is a real parameter. We consider p odd integer, p = 1, 3, 5, ... The case p — 1 is trivial, but in general the new 
potential has the form 



This potential was already investigated in Sec. II. The first-order equations are 



dip 
dx 



and their solutions can be written as 



ip±(x) = ±tanh^(a;) 



as it can be shown directly. 

Another example is given by the potential 



(184) 

(185) 
(186) 

(187) 



We can use the deformation (p tanh(x) to obtain 



"^(x) = ^sech (x) 



(188) 
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Both models were investigated in Sec. II. In the second case, the solutions x±{x) = ±arcsinh(a;) can be deformed to 
give 

'^±W=±^^^ (189) 
Vl + a;^ 

which solve the first model. 

It is interesting to notice that the potential (187) has no minima, and their topological solutions connect inflection 
points. A similar result was presented in Ref. [69], where the authors find topological solutions which connect an 
inflection point to a minimum of the potential there considered. 



V. DEFECTS IN HIGHER SPATIAL DIMENSIONS 



There is an old result [70, 71], which informs that models described by a single real scalar field cannot support 
stable solutions in space dimensions higher than 1. This result is known as the Derrick theorem, but we call it the 
Derrick- Hobart theorem, since it was also obtained by Hobart in Ref. [70]. We can show this result as follows: in 
{D, 1) space-time dimensions, let us suppose that there is a model described by some non-negative potential, in the 
form 



jC=-d^<Pd^cj>-V{cl>) 



We also suppose that (j){x) is static solution, with energy 

E = Eg + Ep = j dx(^Vcj>-Vcj> + V{cj>)^ 



Now, for 4>^{x) = (j){\x) we have 



= j dx Qv0^ • y4>^ + v{(i)^)^ 



X^-'^Eg + X-^E^ 



Thus, we can write 



dE^ 



d\ 



A=l 



{2-D)Eg-DEp 



= {2 - D){1 - D)Eg + D{1 + D)Ep 



(190) 



(191) 



(192) 

(193) 
(194) 



and so, for E^ to have a minimum at A = 1, we have to have {2-D)Eg = DEp and {2 — D)(l — D)Eg + D{l + D)Ep > 0, 
and this requires that D = 1. 

In spite of this result, in a recent work we have been able to circumvent the standard situation, to show how to 
construct stable defect solutions for D arbitrary [35] . To make this possible, we work with the model 



where 



/: = -d^ci>d^(t>-u{x,ci>) 



U{x,<f>)=F{x^x'')V{<f>) 



(195) 



(196) 



and V{(t>) is some standard potential, and F{xij,x^) is a new function, which depends only on the space-time coordinates 
in {D, 1) dimensions. See Ref. [72] for other investigations concerning higher dimensional defect structures. 
In the case of static solutions we can choose an appropriate F, to write the equation of motion in the form 



(197) 
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Thus, for V = (1/2) W| and for radial solutions, ^ = ^{r) should obey 

n_i d / n_irf(/>\ 



We can write 



to change the equation of motion to 



dr 



dr dy 



= w,w,. 



(198) 



(199) 



(200) 



This equation has the form of the equation of motion for static solutions in the case of two-dimensional models. Thus, 
it can be written as first-order equations 



dy 



or better 



dr 



(201) 



(202) 



This result can also be obtained from the energy of static, radially symmetric solutions. If we work with the energy, 
we can write 



E = -flD 
2 



„D-1 



dr 



dr 



1 



r.2D-2 



, dW 
dr— + 
dr 



Jo 



dr 



(203) 



where Q^D = 27r^/^/r(_D/2) stands for the angular factor, which depends on D. We see that the energy is minimized 
to the value Ebps = ^d^W, where 



AW = W[ 



,)] _ w[(k{r ^ 0)] 



(204) 



if the static, radial solution (j){r) obeys one of the two first-order equations (202). 

As we have shown in Ref. [35], other results can also be obtained. For instance, we can show that the radial solutions 
arc stable against radial fluctuations. Also, wc can introduce fcrmions, with the standard Yukawa coupling changed 
to F((/>) = r^~^Wcj,(j,, to include the radial modification we have introduced in the scalar potential. 

On the other hand, in the case D = 2, the transformation which maps the problem into a one- dimensional problem 
has the form y = ± ln(r). which very much reminds us of the Colc-Hopf transformation, used to transform the Burgers 
equation into a linear, exactly solved equation [73, 74] - see also Refs. [75, 76]. 

To illustrate the general situation, let us consider the model 



t/(x,<^) = ilsin2(0) 
2 



(205) 



in D = 2. This is the sine-Gordon model, but here it is modified to incorporate the spatial dependence introduced 
above. 

In this case, the equation of motion which describes static and radial solutions gets to the modified form 



dr^ ' r dr 

It admits scale invariance, and it can be written as 



d'^cj) Ids 1 . , , ,, 
H ;- = — sm(0) cos((/)) 



dr \ dr 



^31 r 37 = ^H4') cos(0) 



(206) 



(207) 
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FIG. 17: Topological soliton in D = 2, which appears in an isotropic uniaxial ferromagnetic system. 



Thus, we take y = ln(r) to obtain 



— = s\n{4>) cos(0) 



This problem was already solved in Sec. II. We can write 



dy 



= ± sm{(j)) 



and the solutions are 



We use the radial coordinate to get 



:{y) = ±2arctan(e=^^) 



(j)^{r) = ±2arctan(r='=^) 



(208) 



(209) 



(210) 



(211) 



These solutions are topological solutions [35]. They can be used to describe solitons in bidimensional, isotropic, 
uniaxial ferromagnetic [77], and they were first described in Ref. [78]. In Fig. 17 we plot one of the solutions, to show 
the topological structure explicitly; see Ref. [25] for related investigations. 



VI. WARPED GEOMETRY WITH AN EXTRA DIMENSION 



We can couple scalar fields with gravity. We get [79, 80] 



-I 



dtdxy 



\9\ [ 



(212) 
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Here we are working in arbitrary (D, 1) space-time dimensions, using dimensionless fields and variables with AnG = 1. 
(/jui/ is the metric tensor and R is the curvature scalar, given by 

R = i?/ (213) 

where R^i, is the Ricci tensor, which can be written as 

-R/ii/ = R^xu^ (214) 
in terms of the Riemann tensor R^jjx' ^ which can be written in terms of the Christoffel symbols, in the form 

with 

r'^.A = Id^^idugXa + dxg.a - OaQux) (216) 



Einstein equation is given by 



Guv — R/iv — 2^9nv — "^Tfj^v (217) 



and G^jy is the Einstein tensor. 

Our interest on scalar fields coupled to gravity is motivated by String Theory. As one knows, superstrings need 
10 space-time dimensions [81, 82]. In Ref. [83], Polchinski has shown how to construct soliton-like brane solutions 
which appear under the factorization of the (9, 1) space-time into AdS^ x S^. In this case, the interesting structure of 
space-time is AdS^ , and we can work with an extra dimension, following the Randall-Sundrum [84] or Karch- Randall 
[85] scenarios, as we have done in Refs. [86-88]. 

We illustrate the investigations working with the AdS5, in (4,1) space-time dimensions, in the Randall-Sundrum 
scenario. The metric for such warped geometry with an extra dimension is given by 

ds^ = e^^ ri^^dx^dx'' + dy'^ (218) 

where 7)^1, = (—1, 1, 1, 1), and y is the spatial coordinate which describes the extra dimension, e^^ is known as the 
warp factor, since it modifies the standard four dimensional Minkowsky metric, described by 77^,^. The matter contents 
of the model is described by real scalar fields. We consider the action 

S = j d^xdy^M Qi? - ^g^,d^'ct>d''ct> - V{ct>)^ (219) 

where we consider a single real scalar field. For more fields, we modify the scalar Lagrangean accordingly. 

We suppose that A = A{y) and that the scalar field is also static, and only depends on the extra coordinate, 
(j) = (j){y); in this case the equations of motion are 

</," + 4A>' = ^ (220) 
A" + \<j)'^ = (221) 

and 

- \^'^ + = (222) 

which is a consistence equation. In the above equations, we are using that A' = dA/dy, etc. 

In the Randall-Sundrum scenario [84], Minkowsky space M4 is embedded in anti-de Sitter space AdS^. In Refs. [89- 
91] the authors introduce very interesting investigations on the subject. For instance, if we consider the potential in 
the form 

m-U^V->^ (223) 
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we can write the first-order equations 



1 dW 



(224) 



and 



A' 



(225) 



which solve the equations of motion; see also Refs. [92, 93] for other details on this. We notice that the first-order 
equation for the scalar field is very much like the first-order equation that we get for models described by a single 
real scalar field in (1, 1) dimensions in flat spacetime; see Sec. II. Thus, we can use those models to describe brane 
solutions in scenarios involving one extra dimension. 

An example was considered in [86], where we have used the model described by 



P 



2p-l 



2p + r 



(2p+l)/p 



(226) 



There we have shown that the parameter p gives rise to an interesting effect, leading to thick brane with internal 
structure, as wc show in Figs. 18, 19 and 20 for p = 1,3,5. We notice that the thickness of the brane increases for 
increasing p, and that the matter energy density opens a gap in its center, suggesting the appearance of internal 
structure, as we show in Ref. [86]. We also notice that p works like the temperature, as it is shown in the model 
considered in Ref. [94]. 





FIG. 18: The kink-like solutions for p = 1, 3, 5. The thickness of the line increases for increasing p. 
In another work, we considered another model, described by [87] 

W{^,x)=cp-l^''-rcl>x^ (227) 

This model gives rise to Bloch branes, in analogy with the Bloch walls that it generates in flat space-time. The 
interested reader should see Refs. [95-97] for other investigations. 



A. Other Scenarios 



We can work with other scenarios, involving two or more extra dimensions. In the case of two extra dimensions, 
the important defect structures are vortices, which are planar topological solutions that appear in models defined in 
(2, 1) space-time dimensions, describing interactions between gauge and charged Higgs fields. 

A model under investigation is described by [98] 



jC = ^Ke^"'^A^d,nA>, + D^ipD^ip - V{\^\'') 



(228) 
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FIG. 19: The warp factor for p = 1, 3, 5. The thickness of the hne increases for increasing p. 
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FIG. 20: The matter energy density for p = 1,3, 5. The thickness of the hne increases for increasing p. 

In this case, the Chern-Simons term presents higher-order derivative, and we are considering D^^ = 9^ — ^ee^luX^'^A'^, 
which describes specific non-minimal coupling between the vector and scalar fields. 
In this case, the equations of motion are 

sx^^u^^{KDA•' + r) = (229) 
dV 

D^D''ip + — = (230) 

where 



J'* = ie[<fiD'^ifi - {Dt'ip)ip] (231) 



is the Noether current. 

The energy density for static configurations can be written in the form 
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where D = V + ien x VAq. This result appears after ehminating the constraint 



B = -- 



2e2 



(233) 



The energy gets the form 



E = 2Tr J rdr [^\{Di ± iD2)ip\^ + 
±2tt j rdrV'^Ao 



2e if A 



where we have set 



The energy can be minimized to 



Ebps = ±27r / rdrV'^Ao 



It is obtained for fields that obey the first-order equations 



B=T|(i-kn 

K 



We choose the fields in the form 



dAp 
dr 

where n = ±1,±2, The boundary conditions are 

g{r ^ 0) 
ao(r ^ 0) 



in9 



Thus, the first-order equations change to 



dg 
dr 
1 dao 
r dr 



g{r)e 

- — [«o(r-) - n] 
er 

0. g{r oo) 1 
n, ao{r oo) 

±-gao 



(234) 
(235) 

(236) 

(237) 
(238) 

(239) 
(240) 



(241) 
(242) 



(243) 

(244) 



The energy is quantized, in the form Ebps = 27r|n|. 

This models contains the Chern-Simons term; thus, J° = —2e'^\Lp\'^B shows that the electric charge is given by 
Q = 2'K{R/e)n. On the other hand, the magnetic fiux is given by $b = ±{'KR/e^){n'^ + \n\). Moreover, if we introduce 
the topological current = e^'^^d,jA\, wc get that Qt = ^{ttR/ e'){n'^ + |n|), which can be identified with the 
magnetic flux. This identification also appears in the standard Chern-Simons model, but there the values are difierent 
[99, 100]: ^b = Qt. = ±(27r/e)n. 

Here, however, we can introduce another topological current, dual to the above one. It has the form = 
e'^^^d'^J^. In the standard Chern-Simons model, the corresponding topological charge vanishes. But in the model 
under investigation, it is given by Qt = ±2n{R/e)n, which is exactly the electric charge of the solution. 

Our interest in vortices is related to the interest in braneworld scenarios with two or more extra dimensions, but 
this is another story. 
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